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The drag and diffusion coefficients of charm and bottom quarks propagating through quark gluon 
plasma (QGP) have been evaluated for conditions relevant to nuclear collisions at Large Hadron 
Collider (LHC). The dead cone and Landau-Pomeronchuk-Migdal (LPM) effects on radiative energy 
loss of heavy quarks have been considered. Both radiative and collisional processes of energy loss are 
included in the effective drag and diffusion coefficients. With these effective transport coefficients 
we solve the Fokker Plank (FP) equation for the heavy quarks executing Brownian motion in the 
QGP. The solution of the FP equation has been used to evaluate the nuclear suppression factor, 
Raa for the non-photonic single electron spectra resulting from the semi-leptonic decays of hadrons 
containing charm and bottom quarks. The effects of mass on -Raa has also been highlighted. 



PACS numbers: 12.38.Mh,25.75.-q,24.85.+p,25.75.Nq 



I. INTRODUCTION 



Energy dissipation of heavy quarks in QCD mat- 
ter is considered as one of the most promising probe 
for the quark gluon plasma (QGP) diagnostics. The 
energy loss of energetic heavy quarks (Q) while prop- 
agating through the QGP medium is manifested in 
the suppression of heavy flavoured hadrons at high 
transverse momentum (pt)- The depletion of high 
Pt hadrons (D and B mesons) produced in Nu- 
cleus + Nucleus collisions with respect to those pro- 
duced in proton + proton (pp) collisions has been 
measured experimentally [l|-|3( through their semi- 
leptonic decays. The two main processes which 
cause this depletion are (i) elastic collisions and (ii) 
the bremsstrahlung or radiative loss due to the in- 
teraction of the heavy quarks with the quarks, anti- 
quarks and gluons in the thermal bath created in 
heavy ion collisions. 

The importance of collisional energy loss in QGP 
diagnostics was discussed first by Bjorken [4]. The 
calculations of elastic loss were performed with im- 
proved techniques [H, @ and its importance were 
highlighted subsequently @, Q m heavy ion colli- 
sions. The collisional energy loss of heavy quarks [9( 
has gained importance recently in view of the mea- 
sured nuclear suppression in the pt spectra of non- 
photonic single electrons. Several ingredients like in- 
clusions of non-perturbative contributions from the 
quasi-hadronic bound state [l(|, 3-body scattering 
effects 
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the dissociation of heavy mesons due 
to its interaction with the partons in the thermal 
medium 12] and employment of running coupling 
constants and realistic Debye mass [01 have been 
proposed to improve the description of the experi- 
mental data. Wicks et al. [H| showed that the in- 
clusion of both elastic and inelastic collisions and the 
path length fluctuation reduces the gap between the 
theoretical and experimental results. 



The energy loss of energetic partons by radia- 
tion is a field of high current interest |15l - [l9| . For 
mass dependence of energy loss due to radiative 
processes Dokshitzer and Kharzeev [5o| argue that 
heavy quarks will lose much less energy than light 
quarks due to dead cone effects 2l|. However, Au- 



renche and Zakharov claim that the radiative pro- 
cess has an anomalous mass dependence [22| due to 
the finite size of the QGP which leads to small dif- 
ference in energy loss between a heavy and a light 
quarks. The mass dependence of the transverse mo- 
mentum spectrum of the radiated gluons from the 
heavy quarks is studied in [23| . They found that the 
medium induced gluon radiation fills up the dead 
cone with a reduced magnitude at large gluon ener- 
gies compared to the radiation from a light quarks. 
For high energy heavy quarks the effects of the dead 
cone, however, reduces because the magnitude of 
the angle forbidden for gluon emission behave as ~ 
heavy quark mass/energy (24[. From the study of 
the mass dependence of the radiative loss it is shown 
in [25| that the very energetic charm (not the bot- 
tom) quarks behave like massless partons. Although 
the authors in (2(| concluded that the suppression of 
radiative loss for heavy quarks is due to dead cone 
effects but it will be fair to state that the issue is not 
settled yet. 

The other mechanism which can affect the radia- 
tive loss is the LPM effect [27j which depends on 
the relative magnitude of two time scales of the sys- 
tem (28j : the formation time (rp) and the mean 
scattering time scale (r c ) of the emitted gluons. If 
Tp > t c then LPM suppression will be effective. 
The LPM effect is built-in in the expression for ra- 
diative energy loss of heavy quarks derived in [23l — 
l25l [29l | . In contrast to those, in the present work 
we will separately introduce the LPM effects in the 
energy loss formula. 

The successes of the relativistic hydrodynamical 
model (see [3fj, [3l| for review) in describing the host 



2 



of experimental results from Relativistic Heavy Ion 
Collider (RHIC) [32[ indicate that the thermaliza- 
tion might have taken place in the system of quarks 
and gluons formed after the nuclear collisions. The 
strong final state interaction of high energy par- 
tons with the QGP i.e. the observed jetquench- 
ing [33|,[34{ and the large elliptic flow (u 2 ) [35|,|36| in 
Au+Au collisions at RHIC indicate the possibility of 
fast equilibration. On the one hand the experimental 
data indicate early thermalization time ~ 0.6 fm/c 
[37| on the other hand the pQCD based calcula- 
tions give a thermalization time ~ 2.5 fm/c [38j| (see 
also [23|). The gap between these two time scales 
suggests that the non-perturbative effects play a cru- 
cial role in achieving thermalization. It has also been 
pointed out that the instabilities [40H43I ] may de- 
rive the system toward faster equilibrium. However, 
the inclusion of such effects also does not reproduce 
small thermalization time. 



The perturbative QCD (pQCD) calculations indi- 
cate that the heavy quark (Q) thermalization time, 
is larger [H, H3| than the light quarks and 
gluons thermalization scale Tj. Gluons may ther- 
malized even before up and down quarks [4(1] . 
In the present work we assume that the QGP is 
formed at time Tj. Therefore, the interaction of the 
non-equilibrated heavy quarks with the equilibrated 
QGP for the time interval r.; < t < can be 
treated within the ambit of the FP equation [13, EH 
i.e. the heavy quark can be thought of executing 
Brownian motion [|3, [H, [4914561] in the heat bath of 
QGP during the said interval of time. Therefore, the 
propagation of a heavy quarks through QGP may 
be treated as the interactions between equilibrium 
and non-equilibrium degrees of freedom. The FP 
equation provide an appropriate framework for such 
studies. Boltzmann transport equation has recently 
been applied to study the depletion of high energy 
gluons due to its elastic and inelastic interactions 
with QGP 



The paper is organized as follows. In the next sec- 
tion the evolution of the momentum distribution of 
heavy quarks in QGP are discussed. In section III 
we address the issues of radiative energy loss with 
dead cone effect. The non-photonic electron spec- 
tra is discussed in section IV. The initial conditions 
and space time evolution have been discussed in sec- 
tion V, section VI contains the discussion on the nu- 
clear suppression and finally section VII is devoted 
to summary and conclusions. 



II. EVOLUTION OF HEAVY QUARK 
MOMENTUM DISTRIBUTIONS 

The Boltzmann transport equation describing a 
non-equilibrium statistical system reads: 



0_ 

dt 



E 



F.V r 



f{x,p,t) 



df 
dt 



(1) 



where p and E denote momentum and energy, V x 
(V p ) are spatial (momentum space) gradient and 
f(x,p,t) is the phase space distribution (in the 
present case / stands for heavy quark distribution). 
The assumption of uniformity in the plasma and ab- 
sence of any external force leads to 



dj_ 

dt 



df 



dt 



(2) 



col. 



The collision term on the right hand side of the above 
equation can be approximated as (see [5(3, HH for 
details): 
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where we have defined the kernels 
Ai — J d 3 kuj(p, k)ki 
Bij = I d 3 kuj(p,k)kikj. 



(3) 



(4) 



for | p |— > 0, Ai — > and — > DSij where 
7 and D stand for drag and diffusion co-efficients 
respectively. The function uj{p,k) is given by 
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where /' is the phase space distribution, in the 
present case it stands for light quarks and gluons, 
v is the relative velocity between the two collision 
partners, a denotes the cross section and g is the 
statistical degeneracy. The co-efficients in the first 
two terms of the expansion in Eq. [3] are compara- 
ble in magnitude because the averaging of ki in- 
volves greater cancellation than the averaging of the 
quadratic term kikj. The higher power of fcj's are 
smaller |47j |. 

With these approximations the Boltzmann equa- 
tion reduces to a non-linear integro-differential equa- 
tion known as Landau kinetic equation: 



df 
dt 
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The nonlinearity is caused due to the appearance of 
/' in Ai and B^ through w{p, k). It arises from the 



3 



simple fact that we are studying a collision process 
which involves two particles - it should, therefore, 
depend on the states of the two participating parti- 
cles in the collision process and hence on the product 
of the two distribution functions. Considerable sim- 
plicity may be achieved by replacing the distribution 
functions of one of the collision partners by their 
equilibrium Fermi-Dirac or Bose-Einstein distribu- 
tions (depending on the statistical nature) in the 
expressions of Ai and B^ . Then Eq. [6] reduces to a 
linear partial differential equation - usually referred 
to as the FP equation describing the interaction of 
a particle which is out of thermal equilibrium with 
the particles in a thermal bath of light quarks, anti- 
quarks and gluons. The quantities Ai and By are 
related to the usual drag and diffusion coefficients 
and we denote them by 7^ and Dij respectively (i.e. 
these quantities can be obtained from the expres- 
sions for Ai and B^ by replacing the distribution 
functions by their thermal counterparts). 

The evolution of the heavy quark momentum dis- 
tribution (/) while propagating through QGP can 
be studied by using the FP equation (see [5(| for 
details) 
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During the propagation through the QGP the heavy 
quarks dissipate energy predominantly by two pro- 
cesses: (i) collisional, e.g. gQ — > gQ, qQ — > qQ and 
qQ — > qQ and (ii) radiative processes, i.e. when the 
heavy quark emits gluons due to its interaction with 
the thermal partons in the plasma. Therefore, the 
drag and diffusion coefficient should include these 
two processes of energy dissipation. 

The elastic collisions of heavy quarks with light 
quarks (q) and gluons (g) i.e.: gQ — > gQ, qQ — > qQ 
and qQ — > qQ have been used to evaluate the 
transport coefficients (j C oii and D co u) due to col- 
lisional process. At LHC energy one can not ignore 
the radiative energy loss, therefore, this should also 
be taken into account through the transport coef- 
ficients. The transport coefficient [58j . q, which is 
related to the energy loss [59j . dE/dx of the propa- 
gating partons in the medium, has been used to cal- 
culate the shear viscosity to entropy density ratio, 
rj/s [H, [60|. The q is closely related to the diffusion 
coefficient D (for detail see [60]). In similar spirit 
we use dE / dx to calculate the drag coefficient of the 
medium and use Einsteins relation, D = TMj to 
obtain the diffusion co-efficient when a heavy quark 
of mass M is propagating through the medium at 
temperature T. The action of drag on the heavy 
quark can be defined through the relation: 



dE 

dx 



\rad = IradP 



(8) 



where "/ ra d denotes the drag-coefficient and p is the 
momentum of the heavy quark. It should be men- 
tioned here that the collisional and the radiative 
processes are not entirely independent, i.e. the colli- 
sional process may influence the radiative one, there- 
fore strictly speaking dE/dx and hence the trans- 
port coefficients for radiative and collisional process 
should not be added to obtain the net energy loss 
or net value of the drag coefficient. However, in the 
absence any rigorous way, we add them up to obtain 
the effective drag co-efficients, 7 // = jrad+lcoii and 
similarly the effective diffusion coefficient: -D c // = 
Dcoii + D la d- This is a good approximation for the 
present work because the radiative loss is large com- 
pared to the collisional loss at LHC. With these ef- 
fective transport coefficients the FP equation reads: 
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where 70// and D c fj contain contributions from 
both the mechanisms (collisional and radiative). In 
evaluating the drag co-efficient we have used temper- 
ature dependent strong coupling, a s from [(31) . The 
Debye mass, ~ g{T)T is also a temperature depen- 
dent quantity used as cut-off to shield the infrared 
divergences arising due to the exchange of massless 
gluons. 



III. ENERGY DISSIPATION PROCESSES 

The matrix element for the radiative process (e.g. 
Q + q — > Q + q + g) can be factorized into an elas- 
tic process + Q + q) and a gluon emission 
(Q — > Q + g). The emitted gluon distribution can 
be written as 162:, : 



dn n 
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dr\d 2 k\ 



kj_(kj_ - q A 



tF 2 



(10) 



where k = (k$, k±, k^) is the four momenta of the 
emitted gluon and q = (qo,q±,q3) is the four mo- 
menta of the exchanged gluon, 77 = l/2ln(ko + 
kz)/(ka — k%) is the rapidity and Ca = 3 is the 
Casimir invariant of the adjoint representation and 
a s = g 2 /47r is the strong coupling constant. 

The effects of quark mass in the gluon radiation is 
taken into account by multiplying the emitted gluon 
distribution from massless quarks by F 2 , containing 



the effects of heavy quark mass. 
211: 



fci 



F is given by 



uj 2 Qi + ki 



(11) 



where 8 = M /E. 

As the energy loss of heavy quark is equal to the 
energy which is taken away by the radiated gluon, 
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we can estimate the energy loss of heavy quark by 
multiplying the interaction rate A and the average 
energy loss per collision e, which is given by the av- 
erage of the probability of radiating a gluon times 
the energy of the gluon. 

The LPM effects has been taken into account by 
including a formation time restriction on the phase 
space of the emitted gluon in which the formation 
time, tf must be smaller than the interaction time, 
t = A -1 . The radiative energy loss of heavy quark 
can be given by 



dE , A -l 

-J— \rad= Ae = T \ 

ax 



where e, the average energy per collision is |63, 6 

dn 



(12) 



e = <n g ko> = / drjd k 



d^Ai k ^-^ F2 

(13) 

where Tp = coshr)/k±. As mentioned before for the 
infrared cut-off k™ m we choose the Debye screening 
mass of gluon. 



(14) 



The maximum transverse momentum of the emitted 
gluon is given by: 



fimax\2 ^ (j m ) ^ 

(k± > _< is > 
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FIG. 1: Variation of effective drag coefficient with tem- 
perature for charm quarks 



Following the procedure of earlier works [5CJ, [65J 
we evaluate the drag and diffusion coefficients for the 
elastic processes. Knowing j ra d from the radiative 
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FIG. 2: Variation of effective diffusion coefficient with 
temperature for charm quarks 
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FIG. 3: Same as Fig. |3]for bottom quarks 



processes as described above we obtain the effective 
drag coefficients and hence effective diffusion coeffi- 
cient through Einstein relation. In Figs.[T]and [2]thc 
variation of effective drag and diffusion coefficients 
with T have been depicted for charm quarks. We 
observe that the contribution of the radiative loss is 
large compared to the collisional or elastic one. The 
difference between the collisional and radiative loss 
increases with temperature - indicating very small 
contribution from the former at large T. Similar 
difference is reflected in the diffusion coefficients as 
we have used Einstein's relation to obtain it from the 
drag coefficients. We observe that at low T and pr 
the contributions from collisional processes is more 
than or comparable to that from radiative processes. 
For the bottom quark we find that the gap between 
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the drag coefficients with radiative and eiastic pro- 
cesses is smailer (compared to charm) at lower tem- 
perature domain. Quantitatively the value of drag 
is smaller for bottom than charm quarks because of 
their larger relaxation time. However, the qualita- 
tive behaviour is similar to charm quarks as shown 
in Fig. [31 The diffusion coefficient of the bottom 
quark is large ( Fig]?]) compared to the charm quark 
because of the large mass of the former introduced 
through the Einstein's relation. 

On obtaining the effective drag and diffusion co- 
efficients next we need to know the initial heavy 
quarks momentum distributions to solve the FP 
equation. The production of charm and bot- 
tom quarks in hadronic collisions is studied exten- 
sively [66| . In the present work the pt distribution 
of charm and bottom quarks in pp collisions have 
been taken from the NLO MNR code (67|- The re- 
sults from the code may be tested by measuring the 
production cross sections of heavy mesons (contain- 
ing c and b quarks) in pp collisions at ^/snjv = 5.5 
TeV. With all these required inputs we solve the 
FP equation by using the Greens function technique 
(see [6^j for details). 



IV. THE NON-PHOTONIC ELECTRON 
SPECTRA 



The FP equation has been solved for the heavy 
quarks with the initial condition mentioned above. 
We convolute the solution with the fragmenta- 
tion functions of the heavy quarks to obtain the 
Pt distribution of the heavy mesons (B and D) 
(dN D,B /qTdqr)- For heavy quark fragmentation we 



use Peterson function [68| given by: 

1 



/(*)« 



>[1 



(16) 



for charm quark e c = 0.05. For bottom quark 
e b = (M c /M b ) 2 e c where M c (M b ) is the charm (bot- 
tom) quark mass. The non-photonic single electron 
spectra originate from the decays of heavy flavoured 
mesons - e.g. D — > Xev or B — > Xev at mid- 
rapidity (y = 0) can be obtained as follows [69H70|: 
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q T dq T 



(17) 



where 

F{pT,qr) = w 



d(pT-qT) 
2p T p T .q T 



5 (p T .q T /M) (18) 



where M is the mass of the heavy mesons (D or B), 
uj = 96(1 - 8m 2 + 8m 6 - m 8 - UmHnm 2 )- 1 M~ 6 
(m = Mx/M) and g(E e ) is given by 



9(Ee) = 



E 2 e (M 2 - M x -2ME e f 
(M - 2E e ) 



(19) 



related to the rest frame spectrum for the decay D 
Xev through the following relation [69] 



1 dT H 

T H dE e 



wg(E e ). 



(20) 



We evaluate the electron spectra from the decays 
of heavy mesons originating from the fragmentation 
of the heavy quarks propagating through the QGP 
formed in heavy ion collisions. Similarly the elec- 
tron spectrum from the p-p collisions can be ob- 
tained from the charm and bottom quark distribu- 
tion which goes as the initial conditions to the solu- 
tion of FP equation. The ratio of these two quanti- 
ties, Raa then gives, 



Raa{pi 



dN <! Au+Au 
d 2 p T dy 



N, 



dN' P+P 



(21) 



coll X d' 2 pTdy 



called the nuclear suppression factor, will be unity in 
the absence of any medium. In Eq. [21] N co u stands 
for the number of nucleon-nucleon interactions in a 
nucleus+nucleus collision. The experimental data 
[H-i] at RHIC energy (^5^=200 GeV) shows sub- 
stantial suppression (Raa < 1) for pr > 2 GeV 
indicating substantial interaction of the plasma par- 
ticles with charm and bottom quarks from which 
electrons are originated through the process: c(b) 
(hadronization) — > D(B)(decay) — > e+X. The loss 
of energy of high momentum heavy quarks propagat- 
ing through the medium created in Au+Au collisions 
causes a depletion of high pr electrons. 
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V. SPACE TIME EVOLUTION 

The system formed in nuclear collisions at rela- 
tivistic energies evolves dynamically from the initial 
to the final state. The time evolution such systems 
may be studied by solving the hydrodynamic equa- 
tions: 



Q^U = Q 



(22) 



with boost invariance along the longitudinal direc- 
tion 0- In the above equation T^" = {e+P)u^u u - 
g^ v P, is the energy momentum tensor for ideal fluid, 
e is the energy density, P is the pressure and is 
the hydrodynamic four velocity. It is expected that 
the central rapidity region of the system formed af- 
ter nuclear collisions at LHC energy is almost net 
baryon free. Therefore, the equation governing the 
conservation of net baryon number need not be con- 
sidered here. The radial co-ordinate de pen dence of 
T have been parametrized as in Ref. [521 ] . Some 
comments on the effects of the radial flow are in 
order here. The radial expansion will increase the 
size of the system and hence decrease the density of 
the medium. Therefore, with radial flow the heavy 
quark will traverse a larger path length in a medium 
of reduced density. These two oppositely competing 
phenomena may have negligible net effects on the 
nuclear suppression(see also [52l|). 

The total amount of energy dissipated by a heavy 
quark in the QGP depends on the path length it 
traverses. Each parton traverse different path length 
which depends on the geometry of the system and 
on the point where its is created. The probability 
that a parton is produced at a point (r, <f>) in the 
plasma depends on the number of binary collisions 
at that point which can be taken as: 

P[rA) = ^ {l ~l? )6{R ~ r) (23) 

where R is the nuclear radius. It should be men- 
tioned here that the expression in Eq. (1231) is an ap- 
proximation for the collisions with zero impact pa- 
rameter. A parton created at (r, 4>) in the transverse 
plane propagate a distance L = \J R? — r 2 sin 2 (j) — 
rcos4> in the medium. In the present work we use 
the following equation for the geometric average of 
the integral involving drag coefficient [J cIt"/(t)]: 



J rdrd(j)P(r, </>) / dr^{r) 
J rdrd(/)P(r, </>) 



r = 



(24) 



where v is the velocity of the propagating partons. 
Similar averaging has been performed for the diffu- 
sion co-efficient. For a static system the temperature 
dependence of the drag and diffusion co-efficients of 
the heavy quarks enter via the thermal distributions 



of light quarks and gluons through which it is prop- 
agating. However, in the present scenario the vari- 
ation of temperature with time is governed by the 
equation of state or velocity of sound of the thermal- 
ized system undergoing hydrodynamic expansion. In 
such a scenario the quantities like T (Eq. [24]) and 
hence Raa becomes sensitive to velocity of sound 
(c s ) in the medium. This will be shown in the next 
section. 



VI. THE NUCLEAR SUPPRESSION 

The pt dependence of Raa is sensitive to the na- 
ture of the initial (prior to the interaction with the 
medium) distribution of heavy quarks (67| . For the 
QGP expected to be formed at the LHC we have 
taken initial temperature, Ti = 700 MeV, initial 
thermalization time Ti — 0.08 fm/c which repro- 
duces the predicted hadron multiplicity dN/dy = 
2100 [zl] through the relation: 



2tt 4 



1 



1 dN 



45C(3)4a e// 7ri% dy 



(25) 



where Ra is the radius of the system, £(3) is the 
Riemann zeta function and a e ff = 7r2 ffe///90 where 
9eff (= 2 x 8 + 7 x 2 x 2 x 3 x N F /8) is the degen- 
eracy of quarks and gluons in QGP, Np=immbev of 
flavours We have taken the value of the transition 
temperature, T c = 170 MeV. 

The value of Raa is plotted against the pr of the 
non-photonic single electron resulting from D de- 
cays in Fig [5] The results show substantial deple- 
tion at large pr indicating large interaction rate of 
the charm quarks with the thermal medium of par- 
tons. The sensitivity of the results on the equation of 
state is also demonstrated in Fig. [5] A softer equa- 
tion of state (lower value of c s ) makes the expan- 
sion of the plasma slower, enabling the propagating 
heavy quarks to spend more time interacting in the 
medium and hence lose more energy before exiting 
from the plasma, which results in less particle pro- 
duction at high px- This is clearly demonstrated in 
Fig. [5] It may be mentioned here that c s increases 
with temperature. Therefore, due the higher initial 
temperature of the QGP formed at LHC the value 
of c s may be larger than that of QGP formed at 
RHIC energies. Keeping this in mind we predict 
the nuclear suppression factors for three values of 
c s = l/\/3 (maximum possible), 1 / \/4 and 1/V5 

(FigED. 

The nuclear suppression for the bottom quarks arc 
displayed in Fig. [5] We observe quantitatively less 
suppression compare to charm quarks. The differ- 
ence between the charm and bottom quarks sup- 
pression are affected chiefly by two factors: (i) for 
different values transport coefficients and (ii) for the 
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different kind of initial px distributions. The bot- 
tom quark has less drag coefficients and has harder 
Pt distributions - both these factors are responsible 
for the smaller suppression of bottom quark. The 
present results on Raa may be compared with those 
obtained in [73| in a different approach. 
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FIG. 5: Nuclear suppression factor, -Raa as a function of 
pr for various equation of state for non-photonic single 
electron resulting form D-mesons decay. 
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FIG. 7: Variation of the ratio of nuclear suppression 
factor, -Raa for charm to bottom quarks as a function of 
Pt- 
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FIG. 8: Comparison of Raa obtained in the present 
work with the experimental data obtained by STAR and 
PHENIX collaboration for y/s^ = 200 GeV. The ex- 
perimental data of STAR and PHENIX collaborations 
are taken from [lj] and @] respectively. 



FIG. 6: Same as Fig. [5]for B mesons. 



In Fig [7] we have plotted the ratio: Raa/Raa as 
a function of px, from where the effect of the mass 
and the role of the nature (soft or hard) of the initial 
Pt distributions can be understood (see also [z3]). 

In Fig. [8] we compare the experimental data ob- 
tained by the STAR [l| and PHENIX @ collabora- 
tions for Au + Au collisions at v / snjv = 200 GeV 
with theoretical results obtain in the present work. 



For the theoretical calculations the value of initial 
and transition temperatures are taken as 400 MeV 
and 170 MeV respectively. The value of initial ther- 
malization time is assumed as 0.2 fm/c. These val- 
ues of initial thermalization time and initial tem- 
perature reproduces the total multiplicity at mid- 
rapidity, dN/dy — 1100. We observe that the data 
can reasonably be reproduced by taking velocity of 
sound c s = l/VE. It should be mentioned here that 
the inclusion of both radiative and elastic losses in 
the effective drag enables us to reduce the gap be- 
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tween the experiment and theory without any en- 
hancement of the pQCD cross section as has been 
done in our previous work (65j . 

So far we have discussed the suppression of the 
non-photonic electron produced in nuclear collisions 
due to the propagation of the heavy quark in the 
partonic medium in the pre-hadronization era. How- 
ever, the suppression of the D mesons in the post 
hadronization era (when both the temperature and 
density are lower than the partonic phase) should in 
principle be also taken into account. The suppres- 
sion of the D mesons in the post hadronization era 
is found to be small [75[ , indicating the fact that the 
hadronic medium (of pions and nucleons) is unable 
to drag the D mesons strongly. 

VII. SUMMARY AND CONCLUSIONS 

We have evaluated the drag and diffusion co- 
efficients containing both the elastic and radiative 
loss for charm and bottom quarks. We found that 
the radiative loss is dominant over its collisional 
counterpart. In the radiative process dead cone and 



LPM effects are taken in to account. With these 
transport coefficients and initial charm and bottom 
Pt distributions from NLO MNR [67jj code we have 
solved the FP equation. The solution of FP equation 
has been used to predict nuclear suppression factors 
to be measured through the semi-leptonic decays of 
heavy mesons (D and B) for LHC conditions. We 
find that the suppression is quite large indicating 
that the heavy quarks undergo substantial interac- 
tions in the QGP medium. The ratio of the suppres- 
sion for D and B quarks has also been evaluated to 
understand the effects of mass on the suppression. 
The same formalism has been applied to study the 
experimental data on non-photonic single electron 
spectra measured by STAR and PHENIX collabora- 
tions at the highest RHIC energy. The data is well 
reproduced without any enhancement of the pQCD 
cross section. 
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